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A comparison is made of the (quenched) light hadron spectrum and of simple matrix elements 
for a hypercubic fermion action (based on a fixed point action) and the clover action, both using fat 
links, at a lattice spacing a ~ 0.18 fm. Renormalization constants for the naive and improved vector 
current and the naive axial current are computed using Ward identities. The renormalization factors 
are very close to unity, and the spectroscopy of light hadrons and the pseudoscalar and vector decay 
constants agree well with simulations at smaller lattice spacings (and with experiment). 
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C<3 '. I. INTRODUCTION 

■ This paper is a continuation of earlier work |l],|2| studying the properties of a particular type of improved action for 
studies of quenched QCD, with fermion-gauge field couplings parameterized by "fat links," and a lattice anomalous 
H 1 magnetic moment "clover" term. One of the actions studied here has fermionic couplings extending over a hypercube. 
J> ' The other action is the fat link clover action. Both actions are 0(a 2 ) improved. Features related to their chiral 
. properties appear to be superior to those of standard discretizations, though this does not correspond to a realization 
of an exact lattice symmetry relation. The goal of this paper is to compare a particular implementation of a fixed point 
(FP) action for quenched QCD in four dimensions to a simpler improved action, via the usual tests of improvement. 
Much more elegant studies have been carried out in two dimensions On the way, the fat link clover action turns 
out to be an attractive alternative. 

The hypercubic action is inspired by the fixed point action program Q , applied to fermion actions for QCD in four 
dimensions J| [|. FP actions are classically perfect, which means that they have the following desirable properties: 

First, their spectrum has no lattice spacing a-dependent corrections of the form a™ for any n. Second, and probably 
more importantly, FP actions satisfy the Ginsparg- Wilson fl0|| remnant chiral symmetry condition, namely that 
Oh, the anticommutator of the propagator with 75 is a local operator. As a result, they suffer no additive quark mass 
renormalization and no multiplicative renormalization of axial currents , and satisfy the index theorem |l2|,[l3| . It 
would clearly be a desirable thing to have a version of such an action which could be used in numerical simulations, 
and this paper describes a candidate action which seems to satisfy all of the properties of an FP action very well, 
though not exactly. 

The particular hypercubic action I will test is one whose free-field limit is constructed by blocking out of the 
continuum. It has the usual clover term, normalized to its tree-level value as described in 0. The gauge connections 
of both the hypercubic action and clover action I study are replaced by APE-blocked [O] links, 
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with (x) projected back onto SU(3) after each step, and (n) — U^{n) the original link variable. I take a — 0.45 
and N = 10 smearing steps, chosen because this fattening produces a very small additive mass renormalization. 

There are two differences between this work and Ref. [Q. The first involves physics: since the completion of that 
work we have discovered that the properties of lattice QCD (specifically the spectrum of real fermion eigenmodes) 
changes as the lattice spacing coarsens, and that above about a max — 0.2 fm the would-be zero modes cannot be 
separated from the doubler modes. Since there is a qualitative change in the underlying dynamics, it does not seem 
to be appropriate to attempt to extrapolate to the continuum physics which involves chiral symmetry from lattice 
spacings greater than a max . The second difference is that in this paper I study FP operators, the analogs of the 
familiar "rotated" operators of the Symanzik program. For the action used in Q, the recursion relation for the FP 
operators is contaminated by a redundant eigenvector. I chose to find a different renormalization group transformation 
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(blocking out of the continuum with a Gaussian blocking function) and a different hypercube action, the "Gaussian" 
hypercubic action, in order to construct the operators. 

The outline of the paper is as follows: After reviewing the properties of FP actions, I discuss the ways that 
approximations to FP actions are imperfect. I summarize the physics features of fat link actions. Then I present 
tests of spectroscopy and of vector and axial current matrix elements, at lattice spacing a = 0.18 fm. I end with some 
conclusions and speculations. 



II. A QUICK REVIEW OF FIXED-POINT ACTION FORMALISM 



To find a FP action for QCD, one begins with a set of fermionic (ip n , ip n ) and gauge field (U^n)) variables defined 
either on a fine lattice or in the continuum. The fine action is defined as 

S = f3S g (U)+^A(U)^j. (2) 

S g is the gauge action, i, j label sites, and A(f7) is the fermion action. Introducing a renormalization group (RG) 
blocking kernel with a pure gauge piece T g and a fermionic piece parameterized by a constant K and a blocking 
function a nbtTl normalized by b 

T - /3T g (U, V) +Kj2&n b ~ ^n(&<„J)(*n t ~ (kw)Vw), (3) 

n b 

one integrates out the fine degrees of freedom to construct an action involving coarse-grained variables ^ nb , $„ b and 
V'' ; ,(»/,.i. 

The renormalization group equation 

e- s ' = ( d^dUe- {T+s] (4) 



has a pure gauge FP at g 2 = (j3 — > oo). In that limit the gauge action dominates the integral; its RG equation is 
given by the same steepest-descent equation as for a pure gauge model 

S FP {V) = mm (S FP (U) + T(U, V)) , (5) 

while the fermions sit in the gauge-field background. The fermion action remains quadratic in the field variables, and 
the transformation of the fermion action is given most easily in terms of the propagator 

(A'(V))-i„, = ±6 nbi< + b 2 a(U)nUA(U))-^a(UZ <n , (6) 

where {U} is the field configuration which minimizes Eqn. [5] for a given {V}. For a blocking factor F, it is useful 
to rescale 6a nb , rl = F z / 2 Sl nbin with the Fourier transform of f2„ b , n for a free theory normalized to f2(g) ~ 1 + 0(q 2 ). 
Then (again in free field theory) the momentum-space FP equation for the propagator becomes 

(A'fa))" 1 = ^ + ^E |n(^)| 2 A(i±^)- (7) 

Another useful quantity is the minimizing field, the value of the fine field (as a function of the coarse field) which 
minimizes the exponential in the Gaussian fermion integral: 

= ^2 A -1 (n - m)ba m ,Fm b A(mb - n b )^(rib) 

m,m b 

(8) 

A related object is the fixed-point field p|Jl5j, a local average of fields which blocks into itself under a RG transfor- 
mation, 

ipFp(n) = w ( n — m)ip(m). (9) 
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The averaging function obeys the recursion relation 

w(Fnt, — m)((m — Fmb) — Xw(rib — nib) (10) 

m 

with A = I/F 3 / 2 carrying the canonical dimension of a fermion field. As the scale of the blocking transformation 
F goes to infinity, the FP field and minimizing field coincide. Correlators involving the FP field have no power-law 
cutoff effects-they are also classically perfect. Thus they play the same role for a FP action as the familiar "rotated 
operator" does in the Symanzik improvement program. 

The easiest way to find the FP field, for an action defined by an RG transformation with a scale factor F, is to 
solve Eqn. [h] self-consistently for w(n). Unfortunately, for the action studied in Ref. with a blocking factor 
F = 2 RGT, A = I/F 3 / 2 is not the leading eigenvalue: The leading eigenvector is w(n) — A(n), with an eigenvalue 
of about 2.56/2 3 / 2 . The phenomenon has been described by P. Kunszt in p6| . In terms of physics this does nothing: 
w(n) — A(n) just contributes a contact term to correlation functions. However, from an engineering point of view is 
a serious problem, because it is quite difficult to determine the FP field by solving the recursive equation for w, and 
pulling out a nonleading eigenvector. 

For a free field action constructed by blocking out of the continuum (F — > 00 or equivalently, m a continuous 
variable in the background of fixed rib), the free FP field is readily constructed using the momentum-space version 
of Eqn. |l Therefore, because I wanted to study FP operators, I decided to construct a new action by following the 
MIT group d]§] an d blocking directly out of the continuum. There is a price associated with this-I can no longer 
solve FP equations for the propagator in a nontrivial gauge field background, as I did in I will follow a hybrid 
approach, of constructing the free field FP action, making it gauge invariant by using connections made of fat links, 
and just tuning the fattening to optimize the chiral properties of the action. The construction of the free action is 
described in the Appendix. 



III. IMPERFECTION 



A. Comparison to the standard tests of imperfection 







An FP action is classically perfect, but an approximation to a FP action is not. It is thus an interesting exercise to 
ask how imperfect an approximate FP action will be. 

One can approximately reconstruct the free field FP action and field by blocking out of the continuum. In Eqn. 
the I = term in the mode sum is the most singular for small q and the next higher order term is the 1/K term. 
Keeping only those contributions, a few lines of algebra yields 

Ma) = 1 -*r-g + m +« M 
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C(n) = — — Kn ' 2 (121 

As Q ~ 1 + q 2 + . . ., the reader will recognize these formulas as rather baroque variations on the usual Symanzik- 
improved action and rotated field - up to an overall normalization constant. Indeed 

1 9m 
Ac ( 1 -D + m --(D 2 -m 2 ))(l -—)) + ... (13) 



C-U + ^f^Xl-fV... (14) 

and so we identify the Wilson r— parameter with 1/K = ra/2. 

The most serious lattice artifacts of a Wilson-like action come from its vertices. In the absence of an explicit 
construction of a FP action in a nontrivial background gauge field configuration, what can one say? First of all, 
a simple expansion of the action in powers of a shows that any hypercubic action made of thin or fat links has 
O(a) contributions to its vertex, just like the Wilson action. These contributions must cancel in perturbation theory 
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between the scalar vertices and an additional anomalous magnetic moment "clover" term, which of course is absent 
in the free theory. 

In principle, the clover term could have any coefficient, but one fact we know is that the spectrum of a FP action 
is classically perfect. One can consider the spectrum of a fcrmion with infinitesimal momentum p in an infinitesimal 
external magnetic field B, E = tjiq + p 2 /2itik — B/ms + • • ■ and constrain vtlb — ran — mo, by hand, if necessary. 
This constraint turns out to be identical to the requirement that the clover term have its usual tree level value. Now 
the calculation just paraphrases the old perturbative result of Heatlie, et. al. ]17j and the result is as expected: the 
spectrum of a hypercubic approximation to a FP action, with the clover term, is improved through 0(a 2 ), and if 
matrix elements are measured using FP operators, they also have no 0(a) discretization errors. 

B. Chiral Properties 

Violations of chiral symmetry in Wilson-like actions are an old story : Write the fermion action as 

A(m) = A + m (15) 

and consider a flavor nonsinglet chiral rotation ip(x) — > exp(ie (:r) ^-75)^(2:), ^(x) — > ip(x) exp(ie a (x)^-j5). The 
Ward identity for this rotation is 



5(0(x 1 , . . .,x n )) 
Se a (x) 



(0(x 1 ,...,x n )d ll jS(x) a ) 
X a 

-(0(xi, . . . , x n )m ip(x)—j 5 ip(x)} 
-<0(xi, . . . ,x n ){^(x)—{ 75 , Ao}V>(z))). 

(16) 

The last term is the lattice artifact. In a generic non-FP action, it will mix with lower-dimensional operators 

^){ 75 , A }^(x) = (1 - Z A )d^ J«(s) + ArmP(xh 5 iP(x) + 0(a) + ... (17) 

to give a nonzero additive renormalization to the quark mass and an overall multiplicative renormalization to the 
axial current. It is also responsible for mixing of opposite parity operators in matrix elements appropriate to, for 
example, Bk- In perturbation theory all of these mixings are due to one- loop graphs and are 0(g 2 ). 

The Ginsparg- Wilson relation eliminates these mixings and renormalizations |l2| (or to be more precise converts 
them into local corrections to the Ward identities), by canceling the propagator connecting the operator to the current 
with the anticommutator against the action term from the anticommutator. 

The ingredient of an approximate FP action most necessary to reproduce this feature of a FP action is a fat 
link. The dominant graphs contributing to Za — 1 and Am contain tadpoles, whose contributions from large g 2 are 
suppressed by the soft qqg vertices of a fat link action fel. The presence of the fat link in the one explicit realization 
of the Ginsparg- Wilson relation, the "overlap action" [n9l 



A = f (1 - (1 - f Ao)/^/|(l - |a | 2 ), (18) 

using a thin link A , can be seen easily by making a hopping parameter expansion in A . 

Finally, can one quantify the size of chiral violations expected for an approximate FP action? Imagine constructing 
a FP action by beginning with the Wilson action and performing a series of factor-of-two RGT's. The Wilson term 
violates the Ginsparg- Wilson relation. It is a dimension-five operator, and so under each blocking step its size decreases 
by a factor of 1/2. After N blocking steps the action will have a range O(N) and the size of the violating operator 
will be 0(exp(— JVln2). This exponential decrease of chiral violations with range of the action seems superficially to 
be the same behavior as is seen with domain wall fermions (where N is the length of the fifth dimension p0[). One 
can easily see this behavior in explicit calculations using free field theory. 
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IV. TUNING FOR CHIRAL BEHAVIOR 



The most visible aspect of bad chiral behavior for Wilson-like fcrmions is the presence of exceptional configurations, 
which arise when the Dirac operator Ao has a real eigenmode at A = — m, a value which happens to coincide with 
the bare mass m at which one attempts to construct a propagator for A + m We have argued in Ref. 0] that an 
action whose gauge connections are fat links has a narrower range of real eigenvalues than a thin link action does. 

1.0 

© 



0.8 



S 

o 



0.6 



C\2 



0.4 



0.2 



0.0 



o + 0.45(10) CI. 2 
o x 0.45(10) H 



© 

<3> 



CD 



J I I L 



X 



Xj 



J I L 



-0.2 



0.0 



0.2 



m 



o 



FIG. 1. m„ and the quark mass vs. bare quark mass for various fermion actions. Data are diamonds and fancy diamonds 
for m% and the quark mass from the PCAC relation respectively, for the fat link clover action of Ref. 2, with the clover term 
boosted by a factor 1.2, and octagons and fancy crosses for the Gaussian hypercubic action. Both actions have fat links with 
a = 0.45, N = 10. 

In Ref. [0 we proposed a method for tuning an action to optimize its chiral behavior. This involved measuring 
the spread of real eigenmodes of the Dirac operator, and varying the amount of fattening and the size of the clover 
term to narrow the spread of real eigenmodes. This is computationally costly and if the amount of fattening is large, 
produces only the unsurprising result, that a clover coefficient close to the tree-level value is the optimum one. 

Here, I take a simpler (though more tasteless) approach: fix the clover coefficient to its tree-level value and vary 
the amount of fattening until the additive mass renormalization, measured naively by extrapolating the squared pion 
mass linearly to zero, appears to be small. This is done at small quark masses, but not at masses which are so small 
that the quenched approximation breaks down due to the presence of unpaired instanton modes p^ |. The results 
are shown in Fig. [l]. Data are diamonds and fancy diamonds for and the quark mass from the PCAC relation 
respectively, for the fat link clover action of Ref. ||, with the clover term boosted by a factor 1.2, and octagons and 
fancy crosses for the Gaussian hypercubic action. Both actions have fat links with a = 0.45, N = 10. A jacknife 
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fit of the four lightest masses of m\ = A(rriQ — m c ) or m q — A(rriQ — m c ) (this is the lattice quark mass from the 
PCAC relation; see Sec. VI. C) gives m c — -0.0226(15) from ml and m c = -0.0248(7) from the quark mass, for the 
hypercubic action, and m c — 0.0023(19) and -0.0002(18) for the fat link clover action. 

To make a controlled comparison of the zero mode spectra of these actions, I show in Fig. |^ the real eigenmode 
spectrum of several of these actions on a family of single instanton configurations. The construction of the instantons 
is described in Ref. p3| . They are SU(2) instantons in singular gauge centered at (L/2 + 1/2, L/2 + 1/2, L/2 + 
1 /2, L/2 + 1/2) in a periodic lattice of size L = 8. The vertical line at p ~ 0.95 marks the smallest radius visible to the 
approximate FP gauge action of Ref. (l) . Diamonds and squares show the usual thin- link Wilson and clover actions 
while the other plotting symbols show fat link actions. As the instantons shrink, the real eigenmodes become more 
positive, and then, when the instanton disappears, they collide with another real eigenmode (a doubler) to produce 
a complex conjugate pair of eigenvalues. 
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P 

FIG. 2. Smallest real eigenmode of the massless Dirac operator vs. instanton size for various fermion actions. Diamonds 
and squares show the usual thin-link Wilson and clover actions while the fancy diamonds show the fat link clover action and 
octagons show the hypercubic action (both with fattening a = 0.45 and N = 10 steps). 

For future reference, recall that an action obeying the Ginsparg- Wilson relation would have a zero eigenmode on a 
topologically nontrivial configuration, but when the instanton "falls through the lattice" the eigenmode should jump 
discontinuously to some positive value. Presumably the more improved an action is, the sharper the rise of the real 
eigenmode as the instanton size shrinks. All fat link actions, both hypercube and clover, satisfy that criterion. Based 
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on our studies with fat link clover actions [EL we believe that this feature, like the small mass renormalization, is 
caused by the fat link, not by features of the fermion action which could be seen in the free- field limit. 

V. NUMERICAL TESTS 
A. Spectroscopy 

I carried out quenched spectroscopy on a set of 80 8 3 x 24 background gauge field configurations with a lattice spacing 
a ~ 0.18 fm using the approximate FP gauge action of Ref. Jj, at /3 = 3.70. As a fiducial test of a simpler action 
I recomputed the spectrum of the fat link clover action of Ref. on the same set of background configurations (in 
order to remove all effects of the gauge fields from comparisons). This action also has a = 0.45, N = 10 APE-smeared 
fat links and its clover term is rescaled by a factor Csw — 1-2. 

The cost of the hypercubic action compared to the clover action is a factor of about 17 per multiplication ip = 
(A + m)x during the iterative construction of the propagators. Both the fat link clover action and hypercubic actions 
appear to require about half the number of iterations as the usual thin link clover action to converge to the same 
residue — (A + ?ti)x| 2 /Ix| 2 j presumably because their high Fourier modes decouple from the gauge fields. Neither 
of these actions developed any exceptional configurations over the studied range (tt/p > 0.54 and 0.64). 

The spectrum analysis is completely standard |jj . Both actions reproduce the results of spectroscopy calculations 
of more conventional actions at smaller lattice spacings. The spectroscopy is tabulated in Tables | and fill 
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0.32 


1.116( 3) 


1.246( 5) 


1.935(11) 


2.040(12) 


0.24 


0.959( 5) 


1.117( 6) 


1.726(10) 


1.862(16) 


0.16 


0.789( 6) 


0.994( 8) 


1.517(13) 


1.735(17) 


0.10 


0.649( 6) 


0.902(13) 


1.370(15) 


1.622(23) 


0.06 


0.541( 7) 


0.851(19) 


1.264(21) 


1.543(29) 


TABLE I. Table of best-fit masses, Gaussian hypercube action, (3 = 3.70 (aT c — 1/4). 


K 


PS 
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A 


0.114 


1.094( 3) 


1.226( 4) 


1.887(10) 


1.984( 9) 


0.116 


0.977( 3) 


1.137( 5) 


1.741 ( 9) 


1.862(12) 


0.118 


0.856( 5) 


1.048( 6) 


1.595(10) 


1.740(16) 


0.120 


0.718( 5) 


0.960( 9) 


1.449(12) 


1.616(25) 


0.122 


0.551( 7) 


0.866(18) 


1.280(19) 


1.555(25) 


0.123 


0.437( 7) 


0.807(28) 


1.190(30) 


1.482(32) 



TABLE II. Table of best-fit masses, fat link clover action, (3 = 3.70 (aT c = 1/4). 
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To show the data graphically, I begin with an Edinburgh plot, Fig. ||, containing j3 — 6.15 staggered fermion data 
pi| and data from the two actions shown here, the fat link clover action and the Gaussian hypercubic action. I 
compare scaling violations in hyperfine splittings by interpolating the data to fixed it/ p mass ratios and plotting the 
N/p mass ratio vs. m p a. I do this at three n/p mass ratios, 0.80, 0.70, and 0.65, in Fig. [|. The bursts are from the 
nonperturbatively improved (thin link) clover action of Refs. |^q] and J26|. The octagons show data from staggered 
fermions at j3 — 6.15. The other plotting symbols show our test actions: They are the hypercubic action of Ref. Q| 
(fancy crosses) and of this work (crosses) and the c sw — 1-2 fat link clover action (squares on (3 — 5.7 Wilson gauge 
actions diamonds on FP gauge configurations). The data indicate that the improved kinetic properties of the 
hypercubic action do not affect hyperfine splittings very strongly. The two outer straight lines in the figure are linear 
fits in (amp) 2 to the staggered and nonperturbatively improved clover data. The inner straight line is a constant at 
the value of the extrapolated a — value of the staggered data. 

The fat link clover action on Wilson gauge configurations has a slightly larger hyperfine splitting than the other 
actions. Could this be an effect of the gauge action on the spectrum? Could it be due to the fact that the Wilson 
gauge action overproduces small instantons compared to the FP gauge action? 
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FIG. 3. Edinburgh plot showing /3 = 6.15 staggered fermions from Ref. 25 (crosses), the fat link clover action (diamonds) 
and the Gaussian hypercube action (octagons). 
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FIG. 4. N/p mass ratio interpolated to ir/p — 0.80 (a), 0.70 (b), 0.65 (c) the nonperturbatively improved (thin-link) 
clover action (bursts), c sw = 1.2 fat link clover action (squares on (3 — 5.7 Wilson gauge actions, diamonds on FP gauge 
configurations), hypercubic action of Ref. 1 (fancy crosses), and Gaussian hypercubic action (crosses). 



Differences appear when one compares the dispersion relations (Fig. ||). The squared speed of light, c 2 = (E(p) 2 — 
m 2 )/p 2 , for p = 27r/8(l, 0, 0), is computed by performing a correlated fit to the two propagators. At larger mass the 
clover action's c 2 falls away from unity while the hypercubic action shows a fairly constant c 2 — 1 ~ 0.025. To achieve 
the same mismatch of c 2 — 1 as the hypercubic action, at the same m/T c , one would have to decrease the lattice 
spacing by about a factor of roughly two (assuming c 2 — 1 ~ 0(a 2 )). 
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FIG. 5. c 2 from Eqn. 30 for the fat link clover action (diamonds) and Gaussian hypercubic action (octagons), for (a) 
pseudoscalar (b) vector and (c) nucleon correlators. Hadron masses are shown in units of the deconfinement temperature 
T c = l/(4o). 



B. Vector Currents 



There are many definitions of conserved currents. The simplest one is constructed by writing |T^.^2| £/ M (n) 
1 + igAfj,(n) + . . . and denning 



Jfi(n) = -i 



SS 



SA^n) 



(19) 



In these actions I have replaced the thin links by fat links, V fi (n) = 1 + igB^(n) + . . . which in linear approximation 
corresponds to 



(20) 



for some smearing function c^ u {r). It is more convenient to replace Eqn. with a "fat conserved current" 



Jfi(n) = -i 



SS 



(21) 
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Are 9 P J M (n) = and d^J^n) = consistent? Yes. In general, has two parts. The first is a term proportional to 
Sftv, which is parity-even, (8^ v cs(r)A flj (r) with cs(r + A) = cs(r — A)), and a term which is proportional to (1 — 8^), 
(1 - 8pv)cA{r)A v (r) with ca(t + A) = -ca(t - A). The "1" terms in 1 - £ M „ are of the form Y^ r ,v CA(r)(Jv{n + r) — 

J v (n + r — v)) which vanishes if is conserved, and so Jp,(n) = ^ r d(r)J^(n + r) for some function d(r). Thus 
conservation of the usual current implies (all rather trivially) conservation of the associated fat current . 

In what follows, the conserved vector current is defined using Eqn. The hypercubic action has a "conventional" 
normalization S = mift'tp + . . ., in contrast to the "kappa" normalization of the clover action S = ipij) — nipDip. In 
either case, the local current is V'T^V' an< i the improved current is (C^lniC^P)- 

Now recall the classification of definitions of the current |27j : 

(f\J^\i)cont = Z JX (f\jX\i) latt + 0(a) + 0(g 2 a) + ... (22) 

where Zjx — F(m)(l + c\ g 2 + . . .) (introducing the field rescaling of the usual Wilson or clover actions) . If the current 
is conserved (C), its Z = 1. If it is improved (I), the 0(a) term is zero. The conserved current of either action is 
not improved. The "rotated" or FP currents are improved but not conserved. It's conventional to rescale the results 
of the clover action so that F(m = 0) = 1. This amounts to rescaling the clover field t/j — ► ^/2Kip. I will present all 
results that way even though it is unnecessary for this method of finding renormalization factors. 
I measured three kinds of matrix elements: 

x _ (n(T)\n(0)) 

71 (7r(T)|Jo(t)*KO)>' [ ' 



R X,Y = MWM^Kkj (OA) 
(n(T,0)\Mt,<l) Y W0,q)) ' 1 ) 



and the ratio 

R > {g) - (our w (<f))- { ) 



Note that (0|Ji|V)) = iirriy/fy gives the vector meson decay constant /y. Naively speaking, R% = Z v x and 

R ' = ZyY I Zyx . 

Eqn. |2J is measured for < t < T with T = 10 on 8 3 x 24 lattices (20 for the hypercube and 40 for the clover 
action), on a subset of the bare quark mass values, and Eqn. ^5] is a byproduct of spectroscopy. The source at t = 
is a Gaussian shell source centered on the origin; the sink at T is a point sink projected onto zero momentum, and 
momentum was injected at the current. I was able to get a signal for q = (27r/8)(l, 0, 0) in addition to q = 0. Fits 
of R^'ff_ are done by a correlated fit to the two three-point functions, fits of R^'i^ are done by single-elimination 
jacknife and fits of Rf' Y by a correlated fit to the two two-point functions. 

The renormalization factors are determined sequentially. I begin with a forward matrix element R% of the conserved 
current: current conservation demands that this ratio be equal to unity, or Zjc = 1. Confirming that measurement, 
I can use Eqn. |2^ or ^4] (at q = 0) to determine ivy's for the other currents. Self consistency for lattice perturbation 
theory demands that the Z's should be process-independent, but in general that will not be true: the 0(a) corrections 
to Eqn. [2^ will make the ratios process-dependent. The absence of process-dependence in the ratios of Z's is a 
measure of improvement. 

By itself, the ratio R% is contaminated by wrap-around in time. The wrap-around mostly affects the two point 
function and it can be corrected for, by fitting the two-point function to a hyperbolic cosine (~ A(exp(— fit) + 
exp(—fi(N t — t)))) and the three-point function to a forward-going exponential ~ A/Z exp(— fiT). 



Tables HI and [V show the results for measurements of Eq. |25| and |24| for the hypercube action, and the same 
results for the fat link clover action are in Tables and VI. Figs. ^ and illustrate the results. All the fat link 



Z-factors are close to unity. In contrast, Z v ~ 0.78 for the nonperturbatively-improved thin link clover action at 
(3 = 6.0- 6.2 §|. 

Tables III and [y] show that the ratio of the conserved current to either the local or improved currents is quite 



different at q = and q ^ 0. This is an artifact of imperfection of the conserved current. Although the ratio of the 
local current to the improved current is also 0(a), its momentum dependence is much smaller. The conserved current 
also does not respond well to being folded over into (0\Ji\p). In Tables |^ and [vj I show only the q — ratio with 
the improved current: it is quite different from the ratio of forward matrix elements. 
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Fig. shows a comparison of different measurements of Z v l/Z v i from two point and three point functions. The 
scatter of the points for the hypercubic action is small and appears to be mass- independent. In contrast, the Z-factors 
of the fat link clover action two-point and three-point functions are quite different, especially at higher mass, although 
they appear to extrapolate into each other at zero quark mass. One could argue that since the local current has 0(a) 
errors, and the improved current is 0(a 2 ), we are just seeing an 0(a) effect and the improved current is the one to 
use. However, the same argument should apply to the hypercubic action, and there the effects are much smaller. If 
it is an 0(a) effect, one will need more than a factor of two reduction in the lattice spacing (or in moa) to make the 
fat link clover action as consistent as the hypercubic action. 
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FIG. 6. Ratios of three-point functions at p = from (a) the hypercubic action and (b) the fat link clover action. Octagons 
label R C ' L , squares R ' 1 , and diamonds R ! ' L . 
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FIG. 7. R ! ' L from ratios of three point functions (octagons q = 0, diamonds q = (27r/8, 0, 0), fancy diamonds, 
q — (27r/8, 27r/8, 0)) and two-point functions (fancy crosses q — 0, crosses q — (2tt/8, 0, 0)). (a) hypercubic action (b) fat 
link clover action. 



am q 


<b(<7 = 0) 


R^(q = 2w/8) 


<'o(9 = 0) 


^(9 = 277/8) 




R^(#=2ir/8) 


0.32 


0.970(1) 


1.01(2) 


1.074(1) 


1.19(1) 


0.906(2) 


0.87(2) 


0.24 


0.971(1) 


1.11(2) 


1.081(1) 


1.24(1) 


0.899(2) 


0.91(1) 


0.16 


0.971(1) 


1.03(3) 


1.088(1) 


1.24(1) 


0.892(3) 


0.86(1) 


0.10 


0.971(1) 


1.42(8) 


1.094(1) 


1.35(2) 


0.868(3) 


0.99(3) 



TABLE III. Table of R*g , Gaussian hypercube action, f3 = 3.70 (aT c = 1/4). 
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K 


R^(q = 0) 


R^(q = 0) 


R^(q = 2w/8) 


0.116 


0.574(9) 


0.933(1) 


0.938(2) 


0.118 


0.621(13) 


0.930(1) 


0.931(2) 


0.120 


0.669(4) 


0.925(1) 


0.925(3) 


0.122 


0.732(7) 


0.920(1) 


0.921(4) 



TABLE IV. Table of R%' Y , Gaussian hypercube action, = 3.70 (aT c = 1/4). 









<'o(9 = 0) 


R^(q = 2tt/8) 


<.o(? = 0) 


R'^(q = 2tt/8) 


0.116 


1.308(17) 


1.349(4) 


1.136(2) 


1.026(4) 


1.2383(1) 


1.301(1) 


0.118 


1.227(3) 


1.449(3)) 


1.009(2) 


1.254(4) 


1.216(1) 


1.26(1) 


0.120 


1.116(3) 


1.262(12) 


1.003(3) 


1.43(5) 


1.149(2) 


1.15(1) 


0.122 


1.080(6) 




0.995(5) 




1.084(1) 





TABLE V. Table of 7?*' y , fat link clover action, j3 = 3.70 (aT c = 1/4). 



K 


R? 1 (q = 0) 


R^(q = 0) 


2*r/8) 


0.116 


0.732(22) 


0.978(1) 


0.973(1) 


0.118 


0.725(35) 


0.986(1) 


0.982(1) 


0.120 


0.714(18) 


0.993(1) 


0.989(1) 


0.122 


0.722(26) 


0.997(1) 


0.994(1) 



TABLE VI. Table of R*/ , fat link clover action, (3 = 3.70 (aT c = 1/4). 



Finally, I can compare l//y to the observed vector meson decay widths, by multiplying the lattice value (given in 
Table VII or VIII ) by the appropriate Z factor (extracted from the q = R n ' correlator- the first or third column 
of Table [il| or |V|) . The result is shown in Fig. ||a. All the currents bracket the <p meson decay constant quite nicely, 
although there is considerable uncertainty in the fat link clover result due to the different renormalization factors. 
Vector decay constants have also been computed using the thin link clover action, with a nonperturbatively tuned 
clover term, at smaller lattice spacings (Wilson gauge action, (3 = 6.0 and 6.2) in Ref. p^ j. I plot a comparison 
between those values and the results of this simulation (showing only the improved operator, to avoid clutter) in Fig. 
|§|b. There does not seem to be much to be gained for this matrix element by going to the smaller lattice spacings. 



mass 


local 


improved 


0.32 


0.188(2) 


0.177(2) 


0.24 


0.203(3) 


0.191(3) 


0.16 


0.226(3) 


0.210(3) 


0.10 


0.239(4) 


0.222(3) 



TABLE VII. Table of lattice 1/ fy for the hypercubic action. 



K 


local 


improved 


0.116 


0.162(2) 


0.159(2) 


0.118 


0.181(2) 


0.179(2) 


0.120 


0.203(4) 


0.201(3) 


0.122 


0.236(4) 


0.235(4) 



TABLE VIII. Table of lattice l/f v for the fat link clover action. 
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FIG. 8. (a) Vector meson decay constants at a — 0.18 fm from the hypercubic and fat link clover actions. Diamonds and 
squares are the local and improved operator with the hypercubic action while the crosses and fancy crosses are the same 
operators with the fat link clover action, (b) A comparison of 1/fv using improved operators with the hypercubic action 
(crosses), the fat link clover action (squares), and the nonperturbatively tuned thin link clover action at /3 = 6.0 (diamonds) 
and 6.2 (fancy diamonds) from Ref. 26. 



C. The Axial Current 



I have performed a nonpcrturbative determination of the axial vector current renormalization factor following the 



classical approach of Martinelli and Maiani ]29| . This begins with the Ward identity Eqn. 16, with the operator 

O = A h u (y)Vp{z) where V p c (z) = %p{z)^- r ypip{z) and A b v {y) — ip(z) ^-7^75^(2/)- Standard manipulations lead to the 
identity 

2m g ££<p>K(£yo)F;(o,o) = -ir bd J2(K d (y>yo)v p c (o,o)) -if aed J2(^ b Myo)A d p (o,o))- (26) 

n y y y 

(P a (n) is the pseudoscalar current ?/>(n)4j- 75^(71).) To pass to a lattice description, all currents are rescaled J — > ZlJl 
and the quark mass is replaced by the unnormalized quark mass from the lattice PC AC relation, m q — > p = m q Zp jZ a 
with 

_ Ey(doA o (y,yo)C(0,0)) 

computed through a (correlated) fit to two two-point functions. The lattice Ward identity then becomes 

- 2p{PAV) = (VV) J- (AA) . (28) 

(I computed Za for the local axial current; the Ward identity for the improved current involves additional contact 
terms.) As for the vector current, I evaluated the two point functions and the three point function with a source on 
t = and a sink on t = yo = 10. I used 40 hypercubic propagators and 80 fat link clover propagators, taking the local 



Zy from the forward current matrix element (rightmost data column of Tables |III| a nd |V|). I only looked at three 
lighter quark masses. Correlated fits to the appropriate two point functions, Eqns. |27| and a correlated fit extracting 
(0|Ao|P5) = /pstops, produced the results shown in Tables IX and [x|. The two errors on Za are from the jacknife 
and from the variation of p and Zy . I use a^/a = 0.4049(37) from Ref. § to connect the lattice spacing a and string 
tension a . In the last column, I quote a physical number using a nominal yfo — 440 MeV. A linear extrapolation 
to the chiral limit in the quark mass produces the bottom line. The agreement with experiment from the hypercube 
action for the pion seems acceptable. Notice that Za shows mild mass dependence and that, in the chiral limit, it is 
also quite close to unity. The clover action shows smaller mass dependence. In contrast, the nonperturbatively-tuned 
thin link clover action Za from the rotated operator is 0.79 at Wilson gauge action j3 — 6.0 p8[ 
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mass 


afps 


zk 


fp S /Vcr 


Us, MeV 


0.24 


0.170(2) 


0.865(1)(1) 


0.3632(55) 


160(3) 


0.16 


0.156(2) 


0.902(2) (2) 


0.3475(60) 


153(3) 


0.10 


0.147(2) 


0.916(2)(2) 


0.3326(55) 


146(3) 







0.966(3) 




134(6) 



TABLE IX. Table of lattice axial parameters for the hypercube action. 



K 


afps 


7 L 




Us, MeV 


0.118 


0.135(2) 


0.851(33)(4) 


0.283(13) 


124(7) 


0.120 


0.134(2) 


0.870(79) (4) 


0.288(30) 


127(13) 


0.122 


0.130(3) 


0.872(28) (4) 


0.280(12) 


123(5) 







0.89(5) 




122(10) 



TABLE X. Table of lattice axial parameters for the fat link clover action. 
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VI. CONCLUSIONS 



Few of the results shown here are surprising. The hypercube action is designed to show improved chiral and kinetic 
properties, and it does. The fat link clover action has poorer kinetic properties, but in the small mass limit they can 
be compensated for by going to smaller lattice spacing. Both actions have vector and axial current renormalization 
factors which are very close to unity. 

It is an open question to me, whether the hypercube action is a practical improvement over the simpler alternative 
of the fat link clover action. The cost of the action is a factor of 17 per computation of ip = A% compared to the usual 
clover action, with a gain that typically only half as many steps of the inversion algorithm are required to construct 
A -1 . The fat link clover action Q has all the gain with no additional cost. Of course, its kinetic properties are 
not improved, and so one must reduce the lattice spacing to improve the dispersion relation and scaling behavior of 
matrix elements. Heavy quark physics with the fat link clover action will suffer from the same kinetic artifacts as the 
standard clover action. 

A potential use for the hypercube action is to construct an action which realizes the Ginsparg- Wilson exactly. That 
is an active area of research [ |30| , with most attempts involving an iteration using the Wilson action as a zeroth-order 
action. The eigenmodes of a GW action lie on or close to a circle. The thin link Wilson action, or even the clover 
action, seem to be poor choices for a zeroth-order action, because their free-field eigenvalue spectrum does not look 
anything like a circle, and because their real eigenmode spectrum in background instanton configurations does not 
look anything like a step function (recall Fig. |^). The fat link stiffens the real eigenmode spectrum on instantons. A 
hypercube action like the one described here looks like a much better choice, provided of course the cost of its use is 
less than the gain in number of iteration steps. 

Finally, dynamical fermion simulations with fat link actions still appear |j3lf to require smaller levels of fattening 
than are needed to substantially improve the chiral properties of Wilson-type actions. But for quenched calculations 
either action seems to me to be superior to a thin link action. 
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APPENDIX: BLOCKING OUT OF THE CONTINUUM 



A particularly interesting choice of blocking kernel for blocking out of the continuum is the overlapping transfor- 
mation 

n( g )=exp(-icV). (29) 

In my work, blocking out of the continuum was approximated by beginning with free Dirac fermions on a fine lattice, 
making one RG step with a blocking factor F = 8, and blocking onto a coarse lattice of size 5 4 sites. 

Truncating the action to a hypercube and requiring good behavior for the resulting free hypercubic action leads to 
the parameter choice c = 0.2 and K — 2.05. By beginning with a bare mass m/F on the fine scale I constructed a 
renormalized trajectory of free massive fermions. The parameterization and dispersion relation was well behaved out to 
very large mass (am ~ 3.0). The optimum K fell smoothly with increasing m, to a value of 1.0 at m = 3. This is shown 
in Fig. |9|. In the actual simulations, the hypercubic parameters were fit to the forms pj(m), \j(m) — a exp(— bra — cm 2 ) 
and the on-site coupling was fixed using the functional relation between the p's, A's, and bare mass m. If negative 



bare masses are needed, I extrapolated exactly as described in Ref. |L|.The action is listed in Table XI. This is of 
course not a unique (or probably even an optimal) parameterization. 





A(r) 


offset 




10 


-0.0725 exp(-0. 7092m - 0.0149m/ 4 ) 


110 


-0.0319 exp(-l. 007m - 0.0421m 2 ) 


1110 


-0.0156 exp(-l. 123m - 0.0803m 2 ) 
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1111 


-0.0080 exp(-1.180m - 0.1169m 2 ) 




Po(r) 






10 


-0.1819 exp(-0.9475m + 0.0031m 2 ) 


110 


-0.0318 exp(-1.031m - 0.0502m 2 ) 


1110 


-0.00897 cxp(-1.017m - 0.1047m 2 ) 


1111 


-0.00295 exp(-0.9307m - 0.1633m 2 ) 



TABLE XI. Couplings of the hypercubic action used in this work. The free action is parameterized as A = A(r) + i~f^p^(r). 





A(r) 


offset 







0.1047 + 0.3287m - 0.0639m 2 


10 


exp(-3.4070 - 0.6776m - 0.0639m 2 ) 
exp(-4.2303 - 0.9708m - 0.07609m 2 ) 


110 


1110 


exp(-4.9495 - 1.0860m - 0.1147m 2 ) 


1111 


exp(-5.6120 - 1.1437m - 0.1509m 2 ) 




po(r) 


offset 




10 


exp(-2.4746 - 0.9150m - 0.09179m 2 ) 


110 


exp(-4.2383 - 0.9947m - 0.0846m 2 ) 


1110 


exp(-5.5061 - 0.9809m - 0.1388m 2 ) 


1111 


exp(-6.6218 - 0.8943m - 0.1969m 2 ) 



TABLE XII. The hypercubic approximation to the FP field for the hypercubic action used in this work. The field is 
parameterized as £ = A(r) + i7 M p M (r). 
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The free FP field was constructed by summing Eqn. g, in parallel with the FP action. Again, I truncated the FP 
field to a hyper cube and parameterized the "offset" terms exactly as for the hypercubic action. These couplings are 
shown in Table XII . Both the action and FP field are made gauge invariant by averaging the gauge connections over 
all the shortest paths, exactly as in Ref. JlJ 

I show a few properties of the free field action: Fig. \wj shows dispersion relations of the hypercubic action on a large 
lattice, for several values of the bare mass. In Fig. [lljl show the eigenmode spectrum of the hypercubic Gaussian 
action (with only the positive imaginary part of the complex eigenmodes). The eigenmodes of the FP action lie on 
a circle of radius K/2 centered at the point (K/2, 0), and the approximate action tracks the circle closely. Finally, 
in Fig. [l^ I display the scalar part of the free fermion propagator, A(r) -1 = (3{r) + i7 M a M (r). The Ginsparg- Wilson 
relation for the propagator probes 1/2(75, A(r)^ 1 } = 7s/3(r). This is shown for massless fermions on an ll 4 lattice 
with antiperiodic boundary conditions. For the real FP action, f3(r) should be proportional to a delta-function at the 
origin. For the hypercubic action it lies about 1-2 orders of magnitude below the value of the Wilson action at all 
nonzero r. 
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FIG. 9. Optimum K defining the RT of the massive action. 
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the Wilson action (octagons). 



[1] T. DeGrand, Phys. Rev. D58 (1998) 09 4503; |ticp-lat/9802012 

[2] 
[3] 



hep-lat/9807002. 



T. DeGrand, A. Hasenfratz, T. Kovacs, 

For studies of the Schwinger model, see: C. B. Lang and T. K. Pany, Nucl. Phys. B513 (1998) 645; F. Farchioni and 
V. Laliena, talk presented at Lattice '97, Nucl . Phys. B (Proc. Suppl.) 63 (1998) 907; F. Farchioni, I. Hip, C. Lang, and 
M. Wohlgennant, talk presented at Lattice 98, |hep-lat /9809049|: jhcp-lat /9812018| F. Farchioni, I. Hip, C. B. Lang, Phys. 
Lett. B443 (1998) 214; W. Bietenholz and I. Hip, [hep-lat/99020i9( 
P. Hasenfratz and F. Niedermayer, Nucl. Phys. B414 (1994) 785. 
U. J. Wiese, Phys. Lett. B315 (1993) 417. 
W. Bietenholz and U. J. Wiese, Nucl. Phys. B464 (1996) 319. 

T. DeGrand, A. Hasenfratz, P. Hasenfratz, P. Kunszt, F. Niedermayer, Nucl. Phys. B (Proc. Suppl.) 53, (1997), 942. 
W. Bietenholz, et al, Nucl. Phys. B (Proc. Suppl.) 53 (1997) 921. 

K. Originos, et al., contribution to Lattice '97, Nucl. Phys. B (Proc. Suppl.) 63 (1998) 904. 
P. Ginsparg and K. Wilson, Phys. Rev. D25 (1982) 2649. 
P. Hasenfratz, Nucl. Phys. B525 (1998) 401. 

P. Hasenfratz, V. Laliena and F. Niedermayer, Phys. Lett. B427 (1998) 125. 
M. Liischer, Phys. Lett. B428 (1998) 342. 

M. Falcioni, M. Paciello, G. Parisi, B. Taglienti, Nucl. Phys. B251[FS13], 624 (1985). M. Albanese, et. al. Phys. Lett. B192, 



21 



pi; 



163 (1987). 

T. DeGran d.A. Hasenfratz P. Hasenfratz, F. Niedermayer, Nucl. Phys. B454 (1995) 587; Nucl. Phys. B454 (1995) 615. 



hep-lat/9706019 Nucl. Phys. B516 (1998) 402. 



P. Kunszt, 

G. Heatlie et al, Nucl. Phys. B352 (1991) 266. 

M. Bochicchio, et. al., Nucl. Phys B262 (1985) 331; L. Karsten and J. Smit, Nucl. Phys. B183 (1981) 103. 

H. Neuberger, Phys. Lett. B417 (1 998) 141. 

Cf. T. Blum, review at Lattice 98, ^iep-lat/9810017 . 

W. Bardeen, A. Duncan, E. Eichten, G. Hockney and H. Thacker, Phys. R ev. D57 (1998) 1633; W. Bardeen, A. Duncan, 
E. Eichten and H. Thacker, Phys. Rev. D 57 (1998) 3890; |hep-lat/9606002 



Cf. R. Mawhinney et al, |hep-lat/9811026] , 
T. DeGrand, A. Hasenfratz, T. Kovacs, Nucl. Phys. B505 (1997) 417. 
Bernard, et al, Nucl. Phys. B(Proc. Suppl.) 53 (1997) 212. 
Gockeler, et. al., Phys. Rev. D57 (1998) 5562. 

G. Edwards, U. M. Heller and T. R. Klassen, Phys. Rev. Lett 80 (1998) 3448; |hep-lat/9711052 
Martinelli, C. Sachrajda, A. Vladikas, Nucl. Phys. B358 (1991) 212. 
Liischer, S. Sint, R. Sommer and H. Wittig, Nucl. Phys. B491 (1997) 344. 
Maiani and G. Martinelli, Phys . Lett. B178 (1986) 285. 



Cf. H. Neuberger |hep-lat /980301 U R. Edwards, U M. Heller, and R . Narayanan, [hep-lat/9807017| , P. Hernandez, K. Jansen 

and M. Liischer, bep-lat/9808010| , and A. Borici, |hep-lat/9810064|. 

Cf. K. Orginos and D. Toussaint, talk presented at Lattice 98, |icp-lat/980914b . 



22 



